April 6, 2010 



Perturbation analysis of trapped-particle dynamics in axisymmetric dipole geometry 

F.-X. Duthoit, A. J. Brizard*, Y. Peysson, and J. Decker 
CEA, IRFM, F-13108, Samt-Paul-lez- Durance, France 

The perturbation analysis of the bounce action-angle coordinates (J, i^) for charged particles 
trapped in an axisymmetric dipole magnetic field is presented. First, the lowest-order bounce action- 
angle coordinates are derived for deeply-trapped particles in the harmonic-oscillator approximation. 
Next, the Lie-transform perturbation method is used to derive higher-order anharmonic action-angle 
corrections. Explicit expressions (with anharmonic corrections) for the canonical parallel coordinates 
s( J, and p|| (J, () are presented, which satisfy the canonical identity {s, p\\}{J, C) = 1- Lastly, ana- 
lytical expressions for the bounce and drift frequencies (which include anharmonic corrections) yield 
excellent agreement with exact numerical results. 

I. INTRODUCTION 

The dynamics of magnetically-confined charged particles exhibits three orbital time scales [l| associated with the 
fast gyromotion, the intermediate bounce motion, and the slow drift motion around, along, and across magnetic-field 
lines, respectively. To lowest order, the fast gyromotion involves circular motion of a charged particle (with mass m 
and charge q) in the plane transverse to a magnetic field line. The radius of gyration (or gyroradius) pg = f j_/|a;g| 
is defined as the ratio of the magnitude of the particle's perpendicular velocity v± = |v_l| and the gyrofrequency 
cjg = qB/mc. When the gyroradius pg is much smaller than the magnetic nonuniformity length-scale Lb (i-e., 
e = pg/Ls ^1), the gyromotion is nearly circular and the action variable pB/ujg associated with gyromotion is an 
adiabatic invariant (here p, = •mv\/2B denotes the magnetic moment). 

The slower bounce and drift motions of a guiding-center are described by the guiding-center Hamiltonian 

. ii + (1) 

and the guiding-center phase-space Lagrangian 

Tg - (^A + p||b) . dX - i/g 6t, (2) 



where A denotes the vector potential generating the magnetic field B = V X A, p|| = mu|| denotes the guiding-center's 
parallel momentum, and the term pB plays the role of a spatial potential in which the guiding-center moves. We note 
that the ignorable gyromotion action-angle pair €{mc/q)pAC,g has been omitted in the guiding-center phase-space 
Lagrangian ^ since it plays no dynamical role in bounce-center dynamics. 

Guiding-center dynamics takes place in a four-dimensional phase space with the magnetic moment p treated 
as an invariant and the gyroangle ^g treated as an ignorable angle. The magnetic moment p and the total energy 
E = p'^^/2'm+pB are thus constant parameters for a given guiding-center orbit. By using the parallel spatial coordinate 

s along a magnetic-field line labeled by the potentials {a, /3) = {y^,y^), for which B = Va X V/? and b = dX./ds, the 
guiding-center Lagrangian ([2]) becomes 

Pg = |; "d/? + Pii (ds + b,dy') - Hg di, (3) 

where A = aV/S and b = Vs -I- biVy'' have been substituted into Eq. ([2|), and summation over repeated indices is 
henceforth implied (unless stated otherwise). The non- vanishing covariant components bi = b • d'K./dy^ guarantee 
that, using b • Vi/* = 0, the magnetic curvature b • Vb = (b • Vbi)Vy^ is non- vanishing. In addition, we note that 
the covariant components may also be expressed as b-i ~ — Vs • V?/'/|Vy*p (no summation), and thus the magnetic 
coordinates (a,/?, s) are in general non-orthogonal (if at least one component bi is non- vanishing and we assumed 
Va-V/3EE0). 
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The present work draws its motivation from some practical aspects of the bounce-center phase-space transformation 
used to obtain the bounce-center phase-space Lagrangian |3H5|] 

fb = — a(ij3 + JdC - Hb dt, (4) 

ec 

where the coordinates {a,f3) denote the bounce-center (e.g., banana-center in tokamak geometry) coordinates, while 
the bounce-center action-angle coordinates (J, C) form a canonical pair (where the bounce angle C is ignorable and the 
bounce action J is an adiabatic invariant). In Eq. the bounce-center Hamiltonian Hb(a, l3,t; J, fj.) is a function 
of the bounce-center coordinates (a, /3) and depends parametrically on the bounce-action J and the guiding-center 
magnetic moment /i. Although the bounce-center phase-space transformation [3] follows the same Lie-transform 
perturbation approach Q used to derive the guiding-center phase-space Lagrangian ^ , there are differences that are 
investigated in the present work. 

The remainder of this paper is organized as follows. In Sec. |lll we review the theory of bounce-center dynamics in 
general magnetic geometry. By using magnetic coordinates (a, /?, s), the lowest-order guiding-center parallel dynamics 
is clearly separated from the slower guiding-center drift motion. From the parallel dynamics, the lowest-order bounce 
action- angle coordinates (J, Q are introduced and the breaking of the lowest-order bounce action by the slow guiding- 
center drift motion is discussed. Next, the general formulation of the bounce-center phase-space transformation is 
presented and the crucial role played by the canonical parallel-dynamics relation (s, p|| ) — > ( J, C) is discussed. While the 
explicit proof of this canonical relation is not known in general magnetic geometry, it is investigated in axisymmetric 
magnetic dipole geometry in Sec. IIIIl In this simple magnetic geometry, the lowest-order bounce action is in fact 
an exact guiding-center invariant (which partially negates the need for the bounce-center transformation). We are 
thus able to explore, in Sec. llVi the canonical parallel-dynamics relation for trapped particles, first, in the deeply- 
trapped (harmonic-oscillator) approximation and, then, by considering anharmonic corrections. Lastly, in Sec. |Vl we 
summarize our work and discuss its application to more complex magnetic geometries (e.g., axisymmetric tokamak 
magnetic geometry). 



II. BOUNCE-CENTER DYNAMICS IN GENERAL MAGNETIC GEOMETRY 

The bounce-center dynamical equations Q describe the adiabatic motion of magnetically-confined charged particles 
in a nonuniform magnetic field, in which the fast bounce-motion time scale in the four-dimensional guiding-center 
motion have been asymptotically decoupled from the two-dimensional slow drift motion across the field lines. 



A. Four-dimensional guiding-center dynamics 



We begin our general discussion of bounce-center dynamics with guiding-center Hamiltonian dynamics in four 
dimensional phase space based on the guiding-center phase-space Lagrangian ([3]). The Euler-Lagrange equations 
obtained from the guiding-center phase-space Lagrangian © are [| 



qJ 



dB . dB\ P\\dbj 



P\\ 




dyi ds ) ^ m ds 



dB _ Pj dh \ 
dy^ m dy^ j ' 

db, dB dbi dB 
ds dy^ dyi ds 



(5) 

(6) 
(7) 



where the two-by-two antisymmetric tensor rf^ has components ry^^ = - 
space Jacobian J , which includes the spatial Jacobian (Va X V/3 • Vs)" 
is 



1 = — ?7^^, and the total guiding-center phase- 
^ — B^^ and the guiding-center Jacobian Bp 



B*. 



cpw - „ - ^ cpii - . / db. , db., 
qB q \dy^ ds 



(8) 



Equation ([S]) describes the slow guiding-center motion across magnetic field lines (in the magnetic- field label-space). 
Equations (jH])-©, on the other hand, describe the fast guiding-center parallel dynamics (at lowest order), with higher- 
order corrections involving magnetic-field curvature and the non-orthogonality of the magnetic coordinates. Note that 
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the guiding-center equations ©-([T]) exactly satisfy the Liouville identity 

which ensures that the four-dimensional guiding-center dynamics conserves the guiding-center phase-space volume. 

The lowest-order guiding-center equations motion are obtained by setting the ordering parameter e = in Eqs. (|6])- 
([5]), so that 2/0 = and 

1 . . 9p\i OB 

So = v\\ and p||o = sq = - ^J■-g^■ (9) 

These equations describe the lowest-order guiding-center parallel dynamics, where the parallel momentum p|| is treated 
in Eq. Q as a function of the magnetic coordinates (a,/3, s) as well as the guiding-center invariants {£,fi): 



\p\\\ia,^3,s■,£,^i) = ^2m [£ - ^iB{a, (3, s)]. (10) 

By conservation of energy £ and magnetic moment fi, a particle following a guiding-center orbit may encounter a 
turning point along a field line where p|| vanishes (at a bounce point s = Sb where sq = 0) and a guiding-center 
becomes trapped between two such turning points s^. 

The lowest-order guiding-center motion Q is periodic on the (s,p||)-plane. The area enclosed by the lowest-order 
periodic orbit can be used to define the bounce action 

j(a,/3;£:,/i) = j Pw^s = \ j ^_ bii K"' ■s; A^) ^^s, (11) 

and the canonically-conjugate bounce angle 

/ N r ds' 

C = TT -(- sgn(p|| Wb / I ,y (12 

Js- \V\\\{S') 

In Eq. (fT2)) . the bounce angle is defined so that ( = ir d.t the turning point and the bounce frequency Wb = 27r/Tb 
is defined in terms of the bounce-period integral 

rb(a,/3;£,M) ^ ^ 2 T ^. (13) 



Vll 



We note that, according to the standard guiding-center time-scale ordering, the ratio a;b/|wg| = 0{e) is small. 



B. Three-dimensional guiding-center dynamics 



We note that if we replace the guiding-center parallel momentum p\\ with the guiding-center energy f , by using 
the relation (I10|) and using the guiding-center energy conservation law £ = 0, we may reduce the description of 
guiding-center dynamics to three coordinates (a,/?, s): 



y 



ri- 



ce 



1 



3_ 

\v\\\ 



9p|| 

w 

ce 



d 
OS 

_ z]_d_ 
q dy^ 



ce 



■ V 



1 



dB 

w 

ce 



— 77 — - 

q dy^ 



Tsi}n\b) 



(14) 
(15) 



where the new total guiding-center phase-space Jacobian is = J'/lwnl, so that the Liouville identity for guiding- 
center dynamics now simply reads 



0. 



The guiding-center motion (jl4p - (ll5l) is thus represented in a three-dimensional phase space (a, /?, s; f , /i), where the 
rapid bounce motion s(C, J) is parameterized by the bounce action-angle coordinates (J, 0, while the slow drift motion 
takes place in the two-dimensional space of magnetic-field labels (a, /3). 
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C. Bounce-action adiabatic invariance 

While the bounce action (jlip is an invariant for the lowest-order parallel guiding-center dynamics (at zeroth order 
in e when the magnetic- field labels a and /3 are frozen), it is not conserved at higher order. Indeed, we find 

dJ ,dJ y' f dBds 

By introducing the bounce-angle averaging procedure 




(where X]± denotes a sum of the sign of v\\) we note that, using Eq. ([T4|) with J = 1 {a± lowest order in e), we find 

where the two-by-two antisymmetric tensor rjij has components ryi2 — 1 — —7721 (i.e., rjij-q^^ = 5^). Equation (fTO)) 
therefore becomes 

which satisfies the standard condition {dJ/dt)h = for adiabatic invariants Q. 

The general theory of adiabatic invariance [7] for charged-particle motion in magnetic fields allows the construction 
a new adiabatic invariant J that is preserved up to second order in e. The resulting asymptotic expansion J = 
Jo -I- eJi + where the lowest-order term Jo = J is given by Eq. ([Tl]) . satisfies the condition for adiabatic 
invariance: 

^ _ dJ _ ^ / di£ dpJi 
~ dt ~ ^\ dt dt 

where we use the expansion for the guiding-center operator d/dt = dt^/ dt + edi/ dt + ■ ■ ■ , with the lowest-order operator 
do/dt = cjb d/d( expressed in terms of the bounce- angle derivative. We can now see that dJ /dt = 0{e^) is satisfied 
at first order in e if the first-order correction Ji is given as [7| 

Ji = - — / ^dC = ^ f (a(/3)b - /3(a)b)dC. (19) 
cjb J dt cuj^ J \ J 

In order to complete the transformation from guiding-center coordinates (a, /?, s( J, C),P|| ( J, C)) to bounce-center co- 
ordinates (a, (3, J, C), we need to find first-order corrections to the remaining bounce-center coordinates (a, /3, Q. 




D. Lie-transform analysis 

The bounce-center transformation [1, Q 

= («,/3,J,C) ^ = (a,AJ,C) (20) 

is a near-identity transformation on the four-dimensional guiding-center phase space: 

^ z" + eGl + ■■■ , (21) 

which is generated by vector fields (Gi, G2, • • • ) where, at n*^-order, the component removes the fast-angle depen- 
dence in the dynamical equation for dz"'/dt. Here, the component G'( = Ji is already given by Eq. ([T9l) . while the 
remaining components {Gi,G^,G'\) are determined by Lie-transform perturbation method. 

As a result of the bounce-center transformation ()21|) . the guiding-center phase-space Lagrangian ([3]) is transformed 
into the new bounce-center phase-space Lagrangian 

Lb = + d5 = e-iPgo + (Fgi - i^iFgo + d^i) + • • • , (22) 
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where the left side is given by Eq. From the first-order (J, C)-coniponents of Eq. ([22]) . we obtain the relations for 
the gauge function Si{a, /?; J, Q: 



ds/dj\ / dSi/dj' 
ds/dC J \ dSi/dC 
while the components G\ ~ {G" , Cf) are determined from the relation 




(23) 



The boimce action-angle components (Gj'jG^), on the other hand, are generally obtained either at higher order in 
the Lie-perturbation analysis or are more directly obtained from the equations of motion. For example, since the 
bounce-center equations of motion dz°' /dt are independent of the bounce angle, we must have 

(tr /dz^K dz" ( dG\ \ 



from which we obtain 



where the bounce-angle-independent part {GDh is computed at higher order. 

Since the bounce-center transformation (1211) is a canonical transformation, its Jacobian is equal to the guiding- 
center Jacobian J'. By using the general relation 



J^J -^Q;-a{jG^) + ••• - (27) 

we easily find that 



dG\ ^ d f ,\ f .\ (dGi , dG\' 



in P\\ - ^-'J' -II - ^ = - Htt + ^ l (28) 



dy' q ' dy' V^n V ^ \ \\ J \ dJ dC T 



where we used Eq. ([15]). 

We now make a few remarks concerning Eqs. ([23|) - (l24l) . First, in order to solve for the components G\, we need to 
solve for the gauge function Si, which is determined from the relations ([^. Second, we can use these relations to 
obtain the identity 

_ d_f ds\ _ d^Si _ _ d_f ds_ 
^ ^ dJ V II dc) ' dCdJ ^ ^ dC V II dJ 

which establishes that the relation (s,p||) — > C) satisfies the canonical identity 

ds dp\\ ds dp\\ , ^ 

^ii>ii ^ ac a7 - a/ " ^- (2^) 

Hence, the solution for 5*1 requires explicit expressions for the parallel-dynamic coordinates as functions of 

bounce action-angle coordinates (J, C)- Since these expressions are generally not known beyond the deeply-trapped 
approximation (in which the bounce motion is represented as a simple- harmonic oscillation between turning points), 
the bounce-center analysis can only proceed forward on a formal basis (beyond the deeply-trapped approximation) by 
taking the canonical identity (1^^ as an axiom of bounce-center transformation theory (3j] . We note that the equivalent 
canonical relation in guiding-center theory is expressed as 



= I - bb, (30) 



where denotes the guiding-center gyroradius vector (which explicitly depends on the gyroangle Cg), the perpendic- 
ular momentum is p_L = m^ldp^/ dC,^, and the perpendicular Poisson bracket { , }x involves derivatives with respect 
to the gyro-action Jg = (mc/g) /x and the gyroangle C,^. The two-dimensional canonical relation pop can be proved 
explicitly because the relations p (JgXg) and p_L(Jg,Cg) are known. 
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III. GUIDING-CENTER MOTION IN AXISYMMETRIC DIPOLE GEOMETRY 

The motivation for the present work is to explore the canonical relation (s,p||) — > (J, C), and the associated canonical 
identity (j29p , in a magnetic geometry simple enough to allow the analysis of bounce motion beyond the deeply-trapped 
approximation. 



A. Axisymmetric Dipole Geometry 

For this purpose, we consider a pure-dipole magnetic field 

M 



B 



73" 



2 cos 9r + sin 9 \ 



(31) 



where spherical coordinates (r, 0, 0) are used and the constant = B^r^ combines the strength of the magnetic field 
Be at the equatorial radial distance [aX — n/2). The strength of the dipole magnetic field at a point (r, 6') in the 
"poloidal" plane (at fixed toroidal angle (f) is 



B{r, 61) = Be -| VT+Scos^i 



while the parallel unit vector is 



b = 



2 cos 9 r + sin I 
Vl + 3 cos2 9 



We note that, in the axisymmetric dipole magnetic geometry represented by Eq. ()31|) . the toroidal angle (j) is an 
ignorable coordinate 

We can also write the magnetic field (pij) in 2-covariant fiux-coordinate representation 



where the poloidal flux ■0 is 



B = V?/; X V0, 



tp (r, 9) = ipc — sin 



Since a magnetic field line lies on a constant-?/' surface ip = tpc, we find that 

r (ro, 9) — Tc sin^ 9 



(32) 



(33) 



(34) 



draws a single magnetic field line on the magnetic surface ip = ipc- This relation yields the distance along a single 
magnetic-field line 



s(re,0) 



sin^ t 



2\/3 



- I Vl + 3cos2 9' sine' d9' = f Vl + 3cos2 9' sm9'd9' 
In fyiTIcos^^- n/ScosS) - VScos^Vl + 3cos' 



(35) 



which is positive for 9 > tt/2 (below the equator), negative for 9 < tt/2 (above the equator), and zero at the equator 
[9 = 7r/2). The length of a magnetic- field line on the constant-?/ surface labeled by re is 



2 + -i= In (2 + V3~ 



= (2.76...) re 



(36) 



The magnetic-field strength on a single line (on a constant-?/) surface labeled by re), on the other hand, is given as 

Bi9) . 5e + ^ (37) 



sm 
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which imphes that all particles are trapped in this pure-dipole field (since B becomes infinite as 6* — or tt). 
The parallel unit vector b associated with Eq. (|32l) is expressed in terms of the coordinates {ip, (j), s) as 

b = Vs + a(V',s)VV', (38) 
where the covariant components {h^,h^^hs) are 

9X 

&0 = b • 9X/9(/> = and bg = 1. We note that, since a = — Vs • W / \W , the magnetic coordinates {ip,(j),s) are 
non-orthogonal, i.e., 

, r f s 2 cos9 \ , „ 

<r,Q) = + 0, 39 

■0 Vl + 3 cos2 d ) 

off the equatorial plane {Q ^ 7r/2). Moreover, the magnetic curvature b • Vb = Vip da/ds is purely perpendicular to 
the magnetic surface. 



B. Guiding-center motion in axisymmetric dipole geometry 

In axisymmetric dipole geometry, the total guiding-center Jacobian ([H]) is simply J — 1 (since b • V X b = 0), with 
= l/|w|||, and the guiding-center drift equations of motion become 



V' = 0, 
= — 

9 



dB d , . 



'as 



(40) 
(41) 



where Eq. (|40p follows from the axisymmetry of the magnetic-dipole field. A trapped-particle orbit is characterized 
by a pitch-angle coordinate ^ = p\\/p which vanishes at the turning points 6^,. Conservation of magnetic moment and 
energy yields the relation 



1 - e (0) 



B{e) B, ' 

and hence, for a fixed value ^o, the turning-point angle Oh is a solution of the relation 



(42) 



(43) 



since ^(6'b) = 0. By using Eqs. (l34l) and ((37)) . the toroidal drift precession angular frequency (|4ip for a guiding-center 
with invariants (ip^S,/!) can be expressed as 



3c £ sin^ 6* (1 + cos2 



qB^rl (1 + 3 cos' 



2Q\2 



2 - 



yi + 3 cos^ e 

sin^ei 



(44) 



where = P\\c/p — ^^1 — B^/E defines the pitch-angle coordinate at the equator. 
The guiding-center parallel equations of motion ©-([T]), on the other hand, become 



P\\ 



dB 
ds ' 



(45) 
(46) 



which immediately imply that the bounce action (jlip is exactly conserved in axisymmetric dipole geometry since 

J(V';f,M) = ^j>P\\ds (47) 

is now explicitly a function of the guiding-center invariants. Lastly, we note that, by Noether's Theorem, the drift 
action in axisymmetric dipole geometry is simply Jd = j c. 



8 



C. Bounce action-angle coordinates 

A trapped-particle orbit can be represented in terms of the canonical action-angle coordinates ([TT|) -([T2 |) . where the 
bounce action is 



and the conjugate bounce angle is 



2tt 



\i\{0,£,c) VTT3cos2^sin6ld6' 



^ reWb r Vl + 3cos2 9' sm9' 



p/mJe- iei(^',eo) 
where ± = sgn(p|| ) and the bounce frequency Wb is defined as 



dJ 



Vl + 3cos'^9sm9 



d£ 2^pJ 

In Eqs. (P5| - ([5lll) . the magnitude of the pitch-angle coordinate 



d9. 



vr+3cos2^ 



sm 



(48) 



(49) 



(50) 



(51) 



is a function of the equatorial pitch-angle coordinate and vanishes at the bounce points 9^^ . When the integral in 
Eq. (|50p is computed numerically as a function of we obtain the following expression for the bounce period Q 



Wb "b "b 



23 
72 



where 



V2re VV27 



■'gel 



(52) 



(53) 



denotes the bounce frequency of deeply-trapped particles (with — and = u/|ajgc|) and the normalized function 
/ (^c) is shown in Fig.[T] Note that, at = 1 (for marginally-trapped particles), the bounce period ((52|) is finite, with 
/(I) = vT8Lc/(27rrc) ~ 1.86. Figure [T] also shows the approximate expression for /(^c) — 1 + 23^^/72 (dash-dotted 
line), which will be recovered from Lie-transform perturbation theory (see Sec. llVSj) . 

Next, using Eq. (|44l) . the bounce-averaged drift frequency is obtained from Eq. (|44l) and is expressed as Q 



'o/b 



(1 ^ COS^ 



d9 



3pc 

2rb rc 7 (l-f 3 cos^ 61)3/2 



2 sm^9 - (l-<eo) Vl + 3 cos2 



sin^e* 



(1 - ^2) Vl + 3 cos2 9 

where the bounce period Tb is defined in Eq. (|52p . The drift period is therefore defined as 

27r 27r , 27r 1 



Td 



Wd 



rid 



.9(6) 



f^d ^ 6^° " 



where 



3 pc w / Po 



2rl 



V\/27 



rib 



(54) 



(55) 



(56) 



denotes the drift frequency of deeply-trapped particles (with = 0) and the normalized function g (^o) is shown 
in Fig. [21 with g{X) = 3/2. Figure [5] also shows the approximate expression g(^c) — 1 + (dash-dotted line), 
which will be recovered from Lie-transform perturbation theory (see Sec. IIVB[) . Using Eq. (|4T|) . we note that the 
bounce-averaged drift frequency (|54p is also defined as 



Wd 



(57) 
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FIG. 1: Plot of the normalized bounce-period function /(^c) as a function of the equatorial pitch-angle coordinate (solid 
line). In the deeply-trapped approximation (dashed line), / (^c) is replaced by unity in Eq. (|52p . The dash-dotted line includes 
the first-order (anharmonic) correction to be discussed in Sec. IIVBI 



1.6 
1.5 
1.4 




'■^0 0.2 0.4 0.6 0.8 1 



FIG. 2: Plot of the normalized drift-period function g (^c) as a function of the equatorial pitch-angle coordinate (solid line). 
In the deeply-trapped approximation (dashed line), g (^c) is replaced by unity in Eq. (|55[). The dash-dotted line includes the 
harmonic and anharmonic corrections discussed in Sec. II V Bl 



For deeply-trapped particles (fo ^ 1), Eq. ([57]) is approximately given as 

from which we recover Eq. (j56p in the deeply-trapped limit (where 5 — )• 1). 

Lastly, when the bounce period l|5^ and the gyroperiod Tg = 27r/|u;g| are compared to the bounce-averaged drift 
period (|55p . we obtain three well-separated orbital time scales Tg ^ Tb <C Td when pc/r^ ^ 1 for all values of (since 
l</,5<2). 



IV. BOUNCE ACTION-ANGLE COORDINATES FOR DEEPLY-TRAPPED PARTICLES 

While the bounce action (|47l) is now an exact invariant for guiding-center motion in an axisymmetric dipole magnetic 
field, we are still committed to investigating the canonical relation (s,P||) — >■ {J,C) in this simple geometry. For this 
purpose, we first consider the lowest-order harmonic trapped-particle motion in the vicinity of the equator (9 ~ 7r/2) 
and, then, consider the first-order anharmonic corrections. 
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FIG. 3: Equatorial-pitch coordinate |^e| versus the normahzed bounce-latitude (Jh/i^- The dashed line represents the linear 
approximation l^c] — 3(7b/\/2- 



We begin by introducing the latitude angle a = 9 — Tr/2 (with \a\ ^ 7i'/2), so that the distance along a field line 
5| becomes 

s = ^sincr + i singer + •••^ = (^a + + ' ' (59) 

i.e., in the vicinity of the equator, |cr| = |s|/''o ^ 1 represents a normalized parallel distance from the equatorial plane. 
In this approximation, the magnetic-field strength ((37|) along a single field line yields the expression 

5,9.2 75 . 4 9 o 63 4 

— = 1 + - sm^ + y sm"* a + ■■■ = 1 + -a^ + —a^ + ■■■ , (60) 

a ^ 3 ^ a + ■■■ . (61) 



while the unit covariant component (j39p is expressed as 



Next, the particle energy reads 



S = §-+,Bjl + la^ + ■■■], (62) 



and hence, close to the equator, the guiding-center moves in the potential (9/2) /i_Bc cr^ of a harmonic oscillator with 
energy £ — fiBg. (We show below that the term fiBc acts as the potential for the bounce-averaged toroidal drift 
motion.) 

Deeply-trapped particles have turning points located close to the equator (i.e., |^c| ^ 1), so that we find 9^ = 
7r/2 ± (Tb, and the equatorial pitch- angle coordinate is 

3ah 11 2 



l^c| = ^^l-^.^ + --j. (63) 

Figure [3] shows the equatorial pitch-angle j^d as a function of the bounce latitude ah obtained from Eq. (|43l) . where 
we clearly see that, for deeply-trapped particles (tTb '''/2), the equatorial pitch-angle |^c| depends linearly on db 
(for up to — 0.6). 



A. Lowest-order harmonic calculations 



We now calculate explicit expressions for the bounce angle-action coordinates p8|) - (l49|) in the deeply-trapped 
approximation. First, in the vicinity of the equator (|ct| ^ 7''/2), Eq. ()5ip yields 

3 

71 



|e|(^,Ce) = — + ... , (64) 
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where terms of first order in tlie ordering parameter 

e,^al^('-^X«l (65) 



have been omitted beyond the "harmonic-oscillator" approximation in Eq. ([64| . The bounce action (|48|) yields the 
lowest-order expression 

From this expression, we use Eq. (|50p to compute the lowest-order bounce frequency 

'^bo = 1 - et — , (67) 



^88 ) V 2 

where we used 



When we combine these two expressions, where we henceforth use Wbo — to lowest order in ct (the first-order 
correction will be used in Sec. IIVB|) . we find that bounce harmonic-oscillator energy 

Jof^b = ^ = ^ = f - /iSe (69) 

is simply the relative energy as measured from the minimum- i? potential /i-Bc- The lowest-order Hamiltonian therefore 
reads 

Ha = fiB, + Jo rib, (70) 

and the lowest-order normalized bounce period V,^ {dHo/dJo)~^ = 1 is shown in Fig. [T] as a dashed line. 

Next, we compute the lowest-order expression for the bounce angle (q from Eq. (|49| . First, using Eq. ((64|) . we 
compute 



Vl + 3cos2 0' sine' , V2 r da' ^/2 
d9 



mO'X.) 3 J_,^ - a'2 3 

and thus we find the lowest-order bounce angle 



TT — arccos 



(71) 



^ |27r- arccos (cr/cTb) (p|| > O) 
1 arccos (cr/o-b) < O) . 

Lastly, the parallel momentum p|| and the parallel field-line coordinate s are thus expressed (to lowest order in et) in 
terms of the action-angle coordinates (Jo, Co) as 

/ \ 1/2 

V\\ = - {^mJoflbj sin Co = P\\o{JoXo), (73) 



and 

1/2 



/ 2Jo y 



cos Co = so(^o,Co), (74) 



where we have used the bounce frequency ([S7)) for deeply-trapped particles. Equations (|73l) - ([71)) represent the standard 
canonical action-angle transformation for a harmonic oscillator and, thus, the lowest-order transformation (p||o, so) — ^ 
(•^o,Co) satisfies the canonical condition (f^: 

r . dso 9p||o dso dp^o 2 /■ , • 2 /■ 1 /^.n 

{^o,P||o}||o = ac^ aj^ - aj; ac^ = ^° + ^° = ^' 

where partial derivatives are evaluated at constant £ and ^o- We note that deeply-trapped particles in axisymmetric 
tokamak geometry [TO] are also described by the harmonic-oscillator representation ([75)l - ([7^ . 
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B. First-order anharmonic corrections 



We now proceed to obtain the first-order (anharmonic) corrections {Ji,(i) for the bounce action-angle coordinates 
{J — Jq + etJi, C = Co + CtCi) through the Hamiltonian Lie-transform perturbation method. By inverting the relation 
(15^)) for s(cr), we obtain 



s 



ft 

3 



which can be inserted in Eq. ((60|) to obtain 

fj.B = iiBe 



39 



The guiding-center Hamiltonian ([T]), written as 

H = HoiJo) 



et Hi (Jo, Co) , 



(76) 



(77) 



(78) 



will therefore have a lowest-order expression Hq{Jq) and a higher-order anharmonic correction iJi(Jo,Co) that will 
break the bounce-action invariance (i.e., Jq = — et dHi/d(o 0). We will henceforth omit contributions of order 
higher than et beyond the harmonic-oscillator dynamics described in Sec. IIV Al 

Using the Lie-transform perturbation method, we perform a transformation to new bounce action-angle coordinates 
(J, C) in order to restore the invariance of the bounce-action, i.e., J = 0{e^). Using the lowest-order bounce frequency 
([57)) . the expression ([77)) becomes 



fiB = fiBc 



m 



1 - et 



et 



13 
12 



where the contribution et comes from the definition of the magnetic moment fi = £ {1 — et f^)/i?c, while the last 
term [with 13/12 — (39/8) (2/9)] represents the anharmonic correction in the magnetic field ([77)) . We therefore have 
the following expression for the first-order perturbed Hamiltonian 



Hi 



m 
"2" 



ills' 



13 
12 



(79) 



where Ki = 9 Jo/(4'Tirg) is defined to satisfy 

dKi 



dJo 



9 Jo 



2 mr^ ' 



which arises from the first-order correction to the simple-harmonic bounce frequency ()67p . Inserting Eqs. ()66|) and 
([7^1 respectively for {Jo) and s (Jq, Co) yields the first-order Hamiltonian 



Hi 



72 

•^0 



9 cos^ Co 



13 
~6 



4 /■ 
cos Co 



■^0 



23 
16 



41 
12 



cos 2 Co 



13 

48 



cos 4Cq 



(80) 



which has bounce-angle-independent and dependent parts. 



1. Anharmonic bounce action-angle transformation 



We extend our description of trapped-particle dynamics in an axisymmetric dipole magnetic field beyond the 
harmonic-oscillator approximation by considering the coordinate transformation Zq = (Jo, Co) ^ Z"' ~ {J,C)- 

Z" = ZS + ^t{Si, Z^}^^^ + ... , (81) 

where the gauge function Si (Jo, Co) is designed to restore the bounce-angle independence to the bounce-center dy- 
namics in the anharmonic approximation, and {,}||o is defined in Eq. ()75p . Note that, unlike the bounce-center 
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transformation (|2ip . the magnetic-label flux coordinates {ip, (j)) are frozen, i.e., only the parallel motion in phase space 
is transformed by Eq. (|5l]) . 

Using the bounce action-angle transformation ()81|) . we build a new bounce- angle independent Hamiltonian 

H{J) = Ho{J) + etHi{J) + ... - ffo(Jo) + et (^ifi(Jo,Co) - {^i, i/o}||^ ( Jo, Co)) + ••• ■ (82) 

From Eq. (|82p. we obtain the following relations at zeroth order 

Ho = fiB, + J17b, (83) 



and first order 



Hi - ((Hi)b + Hi) - r!b (84) 

where we have divided Hi into a bounce-angle-independent part (i?i)b and a bounce-angle-dependent part Hi = 
Hi — (i/i)b. Since Hi must be bounce- angle- independent, we have 

23 J2 
16 mr| 

while the first-order generating function Si{Jo, Co) is defined as 

^i(Jo,Co) = ^/ ^i(Jo,Co) dCo = 2Ami^^r^^ 41sin2Co + y sin4Co) ■ (86) 

From Eq. (j8ip . the new expressions for the bounce action-angle coordinates {J,C)j which include anharmonic correc- 
tions, are therefore 

J = Jo + etlSi, Jo}„ = Jo + = Jo + , (- 41cos2Co + ^cos4Co') , (87) 

I. J ||o OQo 12mS2br^ V 4 y 



C = Co + {^r, Co}„^ = Co - ^ ^ Co + (41sin2Co - f sm4Co) , (88) 



and 



Equations ((87|) - (l88l) can be easily used to show that the bounce action-angle transformation (|8ip is canonical, i.e., 
{C,J}||o = l + 0(e?). 

2. Anharmonic bounce action-angle dynamics 

The new (bounce-angle-independent) bounce-center Hamiltonian now reads 

23 P 

H ^ fiB, + JOb - et -7 — , (89) 

lb mrj 

where the last term denotes the anharmonic correction. The new bounce frequency (with anharmonic correction) is 
therefore evaluated as 

which is exactly equivalent to the direct calculation leading to Eq. ([5^ . Figure [T] shows the normalized bounce 
period /(^e) = f^b/wb ~ [1 — (23/72) (dash-dotted line) with the first-order anharmonic correction, which yields 

excellent agreement with the numerical result (solid line) up to about — 0.6. 
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In order to prove that the new bounce action (|87|) (with Enharmonic correction) is an adiabatic invariant at first 
order in et, we calculate J = { J, H}^ under the perturbed Hamiltonian H — Hq + eti?i: 



J = Jo + Ca 



where we used the lowest-order Hamilton's equations = and 

etJi / . _ 13 



Jo 



6mr, 



2 ^ 41sin2Co + ysin4Coj = - et ^^b ■ 



dJi 



(91) 



3. Anharmonic canonical condition 
Equations ([75)) - (rM)l now give the parallel canonical coordinates (s,p||) (with first-order corrections) as 



1/2 



P\\oiJX) 



s = Sf) + et\ Si, So \ + 
I J ||o 

^ (ILV^' cosC 



+ ^(J-Jo) + ^(C-Co) + ••• = so(J,C) 

OJo OQo 



(92) 



(93) 



We therefore retain the harmonic-oscillator expressions for the parallel coordinates, but with an anharmonic Hamil- 
tonian. Furthermore, one can immediately verify that 



f -I ds dp\\ ds dp\\ .2a, 2 a 



= 1 + O (e?) 



(94) 



i.e., Eqs. ([221)-® are canonical up to O (ef] 



4- Anharmonic- corrected bounce-averaged drift frequency 



Lastly, we can also obtain anharmonic corrections to the bounce-average drift frequency ([57)) as follows. First, we 
express the magnetic-field strength (1501) at constant Ai = tpVc = B^r^ and s as 



^3 



Next, the flux derivative dB/d4> (at constant s) is evaluated as 



.2 



9 



dB 



1 + et 



15 



(95) 



where we have retained the anharmonic first-order correction. Next, by using the expression ()93p for the parallel 
spatial coordinate s{(), the bounce- average drift frequency ([57|) becomes 



3c£ 
qBc rl 

(l - et f 



^t Y^Ce 



1 



15 
18 



(96) 



where we have combined the definition Be — £ {1 ~ et £,c) with the bounce-averaged expression for Eq. (|95p . Figure 
[5] shows the normalized bounce-average drift period g{^e) = ^^d/<^d — [1 — (1/6) C^]"^ (dash-dotted line) with the 
first-order anharmonic correction, which yields excellent agreement with the numerical result (solid line) up to about 

- 0.6. 
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C. Bounce-center coordinates 



The first-order generating vector field for the bounce-center transformation (|20p is calculated using the expressions 
for the deeply-trapped action-angle coordinates. Equations I^E^ and (|^T|) are used, with the lowest-order expressions 
cjj^^ dC, — mds/p\\, to obtain the expression for : 



c 



cp_ f f /dB\ _ dB_ 



q gwb 

where, using Eq. (j95p . we find 



dB\ dB _ 45Jcos2C 

In the lowest-order deeply-trapped approximation, where only the first term in Eq. (|97l) remains, we find 

= 3 — sin2C + • • • = ^ sin2C + • • • , (98) 

where we used Eq. (pT|) for a and Eq. (|56l) for ilj. Here, we easily verify that (Cf )b = 0. 

It was shown in Section |TT] that the fiux coordinate V' is a guiding-center invariant in axisymmetric magnetic- 
dipole geometry, which ensures that = 0. Hence, the bounce-center position X is a purely toroidal shift of the 
bounce- averaged particle position (X)jj, expressed as 

X = (X>^ + eGf(V',J,C) 5X/50. (99) 

We see that the general bounce-center transformation is greatly simplified in the axisymmetric dipole geometry, and 
is consistent with the purely toroidal nature of the magnetic gradient and curvature drifts in this geometry. 



V. SUMMARY 



We now summarize the main results of the present work. The bounce-center transformation that transforms the 
guiding-center phase-space Lagrangian ([3]) into the bounce-center phase-space Lagrangian ^ relies on the existence 
of expressions for the parallel canonical coordinates s{J,() and p\\{J,Q that are required to satisfy the canonical 
identity (|29p . In general magnetic geometry, this canonical identity is taken as an axiom for the bounce-center 
transformation jSj, since the identity is difficult to prove directly beyond the harmonic-oscillator approximation. For 
practical applications in axisymmetric tokamak geometry [lol |. for example, it is a common practice to use expressions 
for s( J, C) and py (J, C) in the deeply-trapped approximation. The primary purpose of the present work was, therefore, 
to investigate the canonical identity (j29p beyond the harmonic-oscillator approximation in an axisymmetric magnetic 
geometry simple enough to allow explicit analytical results. The axisymmetric dipole geometry was chosen for its great 
simplicity (e.g., all confined particles are trapped) as well as its universality in space plasmas [sj and its applications 
in some recent laboratory experiments . 

By using elegant Lie-transform Hamiltonian perturbation methods, we have been able to obtain explicit expres- 
sions (with harmonic and anharmonic contributions) for the parallel canonical coordinates s{J,Q and p||(J, C) in 
axisymmetric dipole geometry, which satisfy the canonical identity By including anharmonic corrections to the 

deeply-trapped particle dynamics, we were also able to obtain explicit expressions for the bounce frequency (|90p and 
the bounce-averaged drift frequency (|96p that yielded excellent agreement (up to — 0.6) with the exact numerical 
results shown in Figs. [T] and [21 respectively, as well as standard numerical fits Q. 

Future work will focus on the bounce-center transformation in axisymmetric tokamak geometry (l]| beyond the 
deeply-trapped approximation, for which the poloidal fiux tp is no longer an invariant, but is replaced with the drift 
invariant -0* = ■0 — py B^, where p|| = p[|/(ma;g) and i?^ denotes the toroidal covariant component of the tokamak 
magnetic field. 
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